Astronomy & Astrophysics manuscript no. newcirc 


February 2, 2008 


(DOI: will be inserted by hand later) 





Flux and spectral variations in the Circinus Galaxy 

S. Bianchi\ G. Matt\ F. Fiore^, A.C. Fabian^, K. Iwasawa^, F. Nicastro"^ 



(N 
O 

o 



> 

cn 
oo 
vn 
a^ 
o 

o 

Oh 
I 

o 



' Dipartimento di Fisica, Universita degli Studi Roma Tre, Via della Vasca Navale 84, 1-00146 Roma, Italy 
^ Osservatorio Astronomico di Roma, Via dell'Osservatorio, 1-00044 Monteporzio Catone Italy 
^ Institute of Astronomy, Madingley Road, Cambridge CB3 OHA, UK 
Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge MA 02138, USA 

Received / Accepted 

Abstract. We report on a dramatic flux (~ 50% increase in the LECS and MFCS band) and spectral variation between two 
BeppoSAX observations of the Circinus Galaxy performed almost three years apart. Through the analysis of all Chandra 
observations available in the archive, including a new DD T observation on May 2001, we show that a high flux state of an 
extremely variable Ultra Luminous X-ray source (CG X-1: 



Bauer et al. 



2001 ), which is within the adopted BeppoSAX source 



extraction region of 2', is the most likely explanation for most of the observed variation. However, the presence of a high flux 
6.7 keV line and the spectral variation of the PDS in the new BeppoSAX data could be partly due to intrinsic variation of the 
nucleus. Comparing the longest Chandra observation and the BeppoSAX one, we find that the long-term flux variability of 
CG X-1 is not accompanied by a significant spectral variability. We also re-analysed the Chandra HFG nuclear spectra and 
report on the presence of a Compton shoulder with a flux of about 20% the line core, in agreement with theoretical expectations 
for Compton-thick matter. 
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1. Introduction 



Galaxies: Seyfert - X-rays: galaxies - X-rays: individuals: CG X-1 - X-rays: 



The Circinus Galaxy hosts one of the closest (3.8 Mpc) and X- 
ray brightest (F2_iokeV ~ 1-5 x 10"" erg cm"^ s"') Seyfert 2. 
The detailed HST images clearly reveal a compact (< 2 pc) 
active nucleus seen through high obscuration, surrounded by 
extended and complex structures ( Wilson et al. 200C ). 

The first X-ray observation was performed during the 
ROSAT All Sky Survey ( [Brinkmann et al.||1994[ ). A reflection 
dominated spectrum was revealed by ASCA, together with a 
prominent neutral iron Ka line and a number of other lines 
from lighter elements (iMatt et al.l 119961). The BeppoSAX ob- 



below ~ 2 keV is probably contaminated by an extended emis- 
sion and/or off-nuclear sources within ASCA and BeppoSAX 
extraction regions. 

This scenario was basically confirmed by Chandra. Two 
different regions are clearly observed: one compact and spa- 
tially unresolved (<15 pc) which corresponds to the nucleus, 
where the reflection spectrum and the iron line is produced; 
the other, extended ov er about 50 pc, where most of the soft 
emission is produced ( Sambruna et al. 200 l4 R)- Furthermore, 
a number of off-nuclear sources were detected, mostly concen- 
trated within 2' of the nucleus ( Smith & Wilson] |2001 ; Bauer 



servation added a precious piece of information, detecting the 
direct nuclear emission above ~ 10 keV, absorbed at lower en- 



et al. 2001; Sambruna et al. 2001a): at least one of them is 



ergies by a Compton-thick (4 x 10 cm ) material, usually 
iden tified with the torus envisaged in unification models (M att 



et al. |1999| ; puainazzi et alj |1999^ 

The line spectru m above ~ 2 keV clearly originates from 



low ionized matter ( Matt et al. 199(: ; Netzer et al 



1998; Sako 



et al. 2000 ). It was shown that it is fully compatible with reflec- 
tion from the inner surface of a mildly ionized torus, the same 
matter likely responsible for the absorption of the nuclear radi- 
ation: this interpretation le ads to an estimate of the inner radius 
of the torus of ~ 0.2 pc (Bianchi et al. 2001). The spectrum 



likely to have contaminated AS CA and SAX observations dur- 
ing its high flux states (see below). 

In this paper we report on a dramatic flux and spectral 
variation detected in a second BeppoSAX observation of the 
Circinus Galaxy performed on January 2001, almost three 
years later than the first one. In order to discriminate be- 
tween a nuclear variation and the contamination of an off- 
nuclear source, we will make extensive use of all Chandra 
observations available in the archive, together with unpub- 
lished Director's Discretionary Time (DDT) data taken on May 



2001. We adopt the names after Bauer et al. (2001) (from now 



Send offprint requests to: 
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Stefano Bianchi, e-mail: 



on BOl) for the two sources discussed in this paper, CG X- 
1 (a2000 = 14M3"'12.^3, ^looo = -65°20'13") and CG X-2 
(azooo = 14''13"'10.^0, ^zooo = -65°20'44"). 
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2. Observations and data reduction 

2. 1. ROSAT 

The Circinus Galaxy was observed once by ROSAT/PSPC in 
1990 within the All-Sky Survey (e.g. [Voges et al]|1996[ ). Later 
on, it was observed five times by ROS ATfHRl (see Table 
Event files were downloaded from the HEASARC archive and 
analysed with HEAsoft 5.1. 

2.2. BeppoSAX 

The Circinus Galaxy was observed for the first time by 



BeppoSAX (Boella et al. 1997) on 1998 Mai'ch 24 and then 



again on 2001 January 7. Data from the first observation are 
discussed by Guainazzi et al. ( 1999| ) and Matt et al. (1999): we 
defer the reader to these papers for details. 

Data from three insti'uments (LECS, MECS and PDS) 
aboard BeppoSAX will be discussed here. The MECS works in 
the 1-10 keV energy band, with an energy resolution of ^ 8% 
and an angular resolution of ^ 0.7 arcmin (FWHM) at 6 keV. 
The LECS covers a larger energy band, extending down to 0. 1 
keV, with characteristics similar to the MECS in the overlap- 
ping band. The PDS is a passively collimated detector with a 
field of view of approximately 1.5x1.5 degrees, working in the 
range 13-200 keV. The effective exposure times of the 2001 
observation were 2.7 x 10^ s (LECS), 5.2 x 10"* s (MECS) and 
3.8 X 10^ s (PDS). 



As discussed in detail in Matt et al. (1999), there is a sec 



ond bright source within 5' from the Circinus Ga l axy i n the 
LECS and MECS images. Following [Matt et al.| ( |1999| ), we 
avoided contamination from this source extracting data for the 
LECS and the MECS in a conservative radius of 2' centered 
on Circinus. A smaller radius implies a large fraction of lost 
photons: for 1.5', a 20% loss of the total counts is measured in 
the MECS. A background was created using blank sky spectra 
from the same region for the detector field of view and then 
subtracted. PDS data have been reduced following the stan- 
dard procedures described in Matt et al. (1997). The resulting 
net count rates are shown in Table [1|. Data were reduced with 
HEAsoft 5.1 and spectra analysed with XSPEC 1 1.1.0. 

2.3. Chandra 

Following the discovery of a large flux and spectral variation 
in the second SAX observation (see Sect. |^, a Chandra DDT 
observation was performed on 2001 May 2. The arcsec spa- 
tial resolution of Chandra allows us to resolve several differ- 
ent sources within the 2' BeppoSAX adopted source extraction 
region. Circinus had been already observed several times by 
Chandra and a systematic study of the variability for all the 
sources present in the field over a period of 1 .5 years is possi- 
ble (Table 111). 

All observations were performed with the Advanced CCD 
Imaging Spectrometer (ACIS-S: Garmire et al., in preparation). 
Only two (they are actually almost contiguous segments) used 
the High-Energy Transmission Grating Spectrometer (HETGS: 
Canizares et al., in preparation). Different configurations of 



subarray windows, resulting in different frame times, were se- 
lected, depending on the aims of the observations. A lower 
frame time reduces significantly the pileup of the nucleus, at 
the expense of a less efficient use of the exposure time and 
consequent lower statistics for all sources. Data were reduced 
with the Chandra Interactive Analysis of Observations soft- 
ware (CIAO 2.2.1), using the Chandra Calibration Database 
(CALDB 2.10). Spectra were analysed with XsPEC 11.1.0, 
while grating data with Sherpa 2.2.1. In the following, errors 
are at the 90% level of confidence for one interesting parameter 
(A^2 ^ 2.71). 



3. Data analysis 

3. 1. Flux 

3.1.1. BeppoSAX 

The net count rates for the new observation are significantly 
higher than fliose of the 1998 one for flie LECS and the MECS 
by a factor of ^ 50%. In contrast, the count rate for the PDS is 
slightly lower in the new observation, even if marginally con- 
sistent with that measured in the old one (Table |l]). 

To avoid possible systematic errors in the effective count 
rates of the new data, we performed a series of tests. First of 
all, we checked that the backgrounds for the two instruments 
were the same in the two observations. Then we looked at the 
lightcurves of the two sets of data in order to understand if the 
calculated mean count rates could be contaminated by short 
periods of anomalous high or low flux. Both tests indicate that 
the flux variation between the two observations is a real effect. 

Another possible source of contamination could be the sec- 
ond source visible in the MECS and LECS images: a rise in 
its flux from the previous observation could affect the Circinus 
flux even within the conservative 2' region adopted. We then 
extracted the spectrum in a radius of 1 ' centered on the second 
source both for the LECS and the MECS. The net count rates 
are indeed higher for both instruments in the 2001 observation, 
but the change is not statistically significant (1.28±0.06xl 0"^ 
against 1.21 + 0.03 x lO^^ for the MECS and 6.7 + 1.7 x 10"^ 
against 5.4+0. 6x 10"^ for the LECS). It seems more reasonable 
that it is indeed the higher flux of Circinus that contaminates 
the second source, and not the contrary. 

The lightcurve of the new data shows a larger variability 
than the old one. In particular, the MECS data in the soft band 
(below 3 keV) seems to vary periodically, as expected if the 
observed excess has to be ascribed to a high-flux state of an 



off-nuclear source, CG X-1, which is strongly periodic (BOl 



see Sect. 4.1). To test this possibility, we folded the MECS 
soft band lightcurve with the period of CG X-1 (27 ks): the 
result, shown in Fig. |l], clearly indicates the presence of such 
a periodicity, with an amplitude between maximum and mini- 
mum which is ~ 50% of the mean value. A similar behaviour is 
found in the LECS data, when folded with the same period: the 
amplitude is of the same order of that observed in the MECS 
data (Fig. |l]). Evidence of the same periodicity is also present 
in the old data, but the variation amplitude is less, indicating 
a contamination from CG X-1 by only a factor of ~ 10%, in 
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LEGS (xlO-2) 


MECS (xlO-') 


PDS 




Flux 2-10 keV (cgs) 


1998-03-13 
2001-01-07 


4.41 ± 0.08 
6.76 ±0.16 


1.050 + 0.009 
1.511 +0.017 


2.01 ± 
1.93 ± 


0.04 
0.04 


1.4 X 10-" 
2.1 X 10-" 



Table 2. The log of the Chandra observations of the Circinus Galaxy, along with 0. 1-10 keV count rates (in units of 10 count 
s ' for nucleus and CG X-1. Flux estimates are not included because the nucleus is largely affected by pileup (see text for details). 



Date 



Instr. 



Frame Time (s) Exp (s) Reference 



Nucleus CG X-1 



2000-01-16 
2000-03-14 
2000-03-14 
2000-06-15 

2000- 06-16 

2001- 05-02 



ACIS-S 
ACIS-S 
ACIS-S 
ACIS-S HETG 
ACIS-S HETG 
ACIS-S 



3.2 
0.4 
3.2 
2.1 
2.1 
3.2 



964 

4974 

23076 

7122 

60222 

4401 



Smith &Wilsor 2001 



Smith &Wilsor 2001 



Smith &Wilsor 2001 



Bauer et al 



Bauer et al 



2001 



2001 



Sambruna et al 



2001a| Jb| 



This paper 



187 
329 
190 
85 ± 1 
89 ±4 
192 ±7 



14 
9 
3 



35 ±5 
339 + 8 
133 ±2 
16 ± 1 
12 ± 1 
67 ±4 




Fig. 1. MECS (top) and LEGS (bottom) soft band (< 3 keV) 
folded lightcurve over a period of 27 ks (7.5 hr) for the January 
2001 observation. The intensities were normalized on the mean 
value. 



agreement with the low flux state in w hich t he source was gen- 
erally observed by Chandra (see Sect. 3.1.2). 



3.1.2. Chandra 

We analysed all the bright sources within 2' of the nucleus (see 
Fig. 1 in BOl) in all Chandra observations (see Table ||). To 
establish if the nucleus has remained constant, it must be noted 



that a direct comparison between observations with different 
frame time is not possible, because of the high pileup the source 
suffers with high frame times. The count rates for the nucleus 
in the three observations with frame time 3.2 s (which span 
a period of 1.5 years) are remarkably constant with one an- 
other. Moreover, the two contiguous observations with frame 
time 2.1 s are also consistent to each other The only observa- 
tion free of pileup is that with frame time 0.4 s (5% according 
to WebPIMMSQ). 

Unfortunately, it is not possible to compare the 2.1s frame 
time observations with the others, because WebPIMMS does 
not support HETGS count rates with frame time other than 3.2 
s. On the other hand, if we use the count rate for the 0.4 s obser- 
vations to calculate the expected count rate with 3.2 s, we get a 
value somewhat lower than the observed one, but the predicted 
pileup is higher than 40%, where WebPIMMS results are prob- 
ably no longer reliable. Therefore, even if, when comparison 
may be done, no evidence for variations is found, it is not pos- 
sible to establish that the nucleus remained constant in all the 
observations, due to the above-mentioned problems. 

All the analysed sources are always much fainter than the 
nucleus but CG X-1, which, on March 2000, reached the count 
rate oftheAGN (Table |). In all the other observations, instead, 
its count rate is much lower than that of the nucleus. The flux 
observed during its peak in March (Fo.i-iokev = 6 x 10"'^ erg 
cm-^ s-') is fully consistent with the LEGS and MECS excess 
in the second SAX observation. 

3.2. Spectral analysis 
3.2.1. BeppoSAX 

Figure ^ shows the data from the new observation superim- 
posed on the best fit for the old one ( |Mattetal.|[T999| ; Guainazzi 



et al. 1999). A huge excess both in the LEGS and MECS is 



evident, together with large residuals around the iron line en- 
ergy and a decrease in the PDS. We subtracted each other the 
two BeppoSAX spectra (PDS excluded) and performed a fit on 



http;//heasarc. gsfc.nasa.gov/Tools/w3piiimis.html 
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Energy (keV) 



Fig. 2. Data from the new BeppoSAX observation superim- 
posed on the best fit for the old one. 



the resulting spectrum. An absorbed powerlaw (F = 2.3;^!!4) 



0.4' 

is a good model (see Table |^), together with a huge (EW=1. 
keV) iron line at 6.67;^q keV. The required column density 
cm"^) is higher than the Galactic one 



(A^H 

(3.0 X lO^i 



-1.7 



10^ 



cm 



-2. 



Bauer et al. 
IT 



for this model is 5.1 x 10" 

-2 c-1 



2001). The 0.1-10 keV flux 



erg cm - s and the line flux is 
^ ph cm"^ s"'. Alternatively, the spectrum can be 
fitted equally well by a multicolour disk blackbody emission 

^ 1 7+0.4 

^ ■ ' -0.3 



7.8+3-2 X 10 



(see Table ||), with an inner-disk temperature of T\ 
keV and A^h = l-O!^^ x lO^^ cm^^. A 6.7 keV line with the 
same flux as for the previous model is still required. 

3.2.2. The Chandra spectrum of CG X-1 

We then analysed the spectrum of CG X-1 in its high flux obser- 
vation, in March 2000 (Table in order to see if the spectral 
shape of the excess in the LEGS and MEGS data is compatible 
with the spectrum of GG X- 1 as observed by Chandra. Because 
of the selected frame time of 0.4 s, this observation is almost 
free of pileup for GG X-1 (5% according to WebPIMMS). We 
get a good fit (x^ - 77/92 d.o.f.) with an absorbed powerlaw, 
A^H = l.l^!" X 10^^ cm~2 an d F = 1 .4+°- j (in good agreement 
with Smith & Wilson 2001). No iron line is required by the 
data: an upper limit to the flux of a 6.67 keV gaussian line is 
4.1 X 10"^ ph cmT^ s"', which is not compatible with that mea- 
sured by SAX. 

A multicolour disk blackbody emission fits the spectrum 
equally well (x^ - 81/92 d.o.f.), with an inner-disk tempera- 
ture of Tin = 2.4+0-4 j^gy ^jjjj ^ g 1+0.9 ^ ^qIX ^^-1 ^ -pj^jg 

kind of spectrum represents emission from an optically thick 
accretion disk around a black h ole, generally a good m odel for 
Ultra-Luminous X-ray sources (Makishima et al. 2000). 

In order to see if the flux variation of this source is associ- 
ated with spectral variations, we extracted the spectra of GG X- 
1 for all Chandra observations. We decided to exclude from 
our analysis the second observation of March 2000, because 
the combination of a high flux state of the source and a frame 
time of 3.2 strongly affects its spectrum with pileup, and the 
shorter one of June 2000, because the source was in the same 
flux state as the contiguous observation, which was longer by a 



factor 10 (see Table ||). In each case, we tried the two models 
described above, a simple powerlaw and a disk black body. The 
results are summarized in Table |[ 

There is indeed evidence of spectral variations, between 
hard and soft states, i n term s either of different photon indexes 
or disk temperatures. BOl found that in the June observation 
there is a possible hardening of the spectrum during the rising 
phase of the period, suggesting that it could be the effect of an 
absorption at soft energies associated with that phase. Because 
of the short exposure times of the other observations, it is not 
possible to confirm this behaviour and ascribe the spectral vari- 
ability to it. The only observation which is long enough is the 
second one performed in March but, because of the high flux 
state and the selected frame time, it is piled up so heavily that 
it lacks significant spectral informations. However, it should be 
noted that the only observation in which the spectrum appears 
soft is the longest one, where we analyse the mean spectral 
shape of all phases of the period. Moreover, this spectrum is 
fully compatible with that of the SAX residuals, as expected 
since the BeppoSAX observation (an exposure time of 51 ks 
distributed over 60 hours) also includes a number of periods of 
GG X-1. 



3.2.3. The Compton shoulder 

We re-analysed the HEG spectrum of the nucleus to look for a 
Gompton shoulder in th e iron Ko- line, which is not reported in 



[Sambruna et al. ( 2001b ), possibly because the response matrix 
included with the version of GIAO available at that time was 
less accurate than the new one. However, this feature is pre- 
dicted on theoretical grounds as a product of Gompton scatter- 
ing of the fluorescent line photons from optically thick matter 



(e.g. |3eorge & Fabiarj 


1991; 


Matt et al 


1991; 


Reynolds et al. 


1994; 


Vlatt 


2002). 


[wasawa et al. 


(1997) 


observed the Gompton 



shoul der in the AS CA spectrum of NGG 1068. Later on, Kaspi 
et al. (2002) detected this feature in the Chandra HETG spec- 
trum of NGG 3783. We expected it to be visible in the 60ks 
Chandra grating observation. 

Indeed, residuals in excess of the red wing of the 6.4 keV 
iron line are clearly visible in the data, if we use a powerlaw 
as the continuum between 5 and 8 keV and put three gaus- 
sian lines at 6.40, 6.67 and 7.04 keV, as required by the data 
(Fig. 1^). We tried to fit this feature with the BOX ID function in 
S her pa, which has a constant value in a given energy range and 
is identically null outside. This is a reas onable appr oximation 
of the theoretical profile (e.g. Fig. 2 in [Vlatl 2002). The best 
fit boundaries for the box are 6.15 keV and 6.40 keV, while the 
flux is 6.2 + 2.0 X 10"^ ph cm"^ g-i^ tjjj^t 25 + 8% that of the 
iron Ka line. If we keep the box boundaries fixed at the val- 
ues expected from theoretical models (see below), the flux of 
the Gompton shoulder is 5.1 + 1.8 x 10"^ ph cm"^ §-1^ jjjat is 
20 + 7% that of the iron Ka line core. 

These values are fully consistent with the theoretical ex- 
pectations. If it originated by reflection from a Gompton-thick 
material (A^h > 10^^ cm"^)^ the first order Gompton shoulder 
is expected to extend from 6.2436 keV to the core of the iron 
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Table 3. Best fit parameters for CG X-1 for the two models adopted in the selected Chandra observations and the BeppoSAX 
spectrum (see text for details). The reported fluxes are relative to the powerlaw model. 



Parameter 


BeppoSAX 


JANOO 


MAROOa 


JUNOOb 


MAYOl 


A'h {cm--) 

r 


2.0!^-^ X 10^2 

9 ri+O'.S 
^■■J_0.4 

15 (35 dof) 


5.6^^^ X 1021 
9 (7 dof) 


i.i+j; 1 X 1022 

1 4+''-2 

77 (92 dof) 


1.3!"-^ X 1022 
2.5 ±0.2 
67 (56 dof) 


8.5!^-^ X 102' 

1 9+04 
1.Z.0 2 

27 (17 dof) 


Nh (cm-2) 
r,„ (keV) 

x' 


1.0+» » X 10^2 
^ - '-0.3 

16 (35 dof) 


5.6!05 X 1021 

2 0+2 2 

9 (7 dof) 


8.1+"^ X 1021 
? 4+04 

81 (92 dof) 


6.8+1"^ X 1021 
1 35+"" 
60 (56 dof) 


6.9!2 » X 1021 
2 9+1-' 
21 (17) 


Flux 0.1-10 keV(cgs) 


5.1 X 10-'2 


5.8 X 10-" 


6x 10^12 


9 X 10^1^' 


1.2 X 10^12 



rti 

1 1 1 1 1 1 1 1 1 1 1 1 1 — 


1 1 1 1 1 1 1 1 1 1 1 1 1 l| 









6.5 

Energy (keV) 



Fig. 3. The Compton shoulder in the Chandra HEG spectrum 
showed as residuals from the best fit including three gaussian 
lines. 



line, with a centroid at 6.31-6.32 keV and a flux of 15 - 20% 



that of the core (Matt 2002) 



As a further test, we also fitted the Compton shoulder with 
XSPEC, adopting a simple gaussian line. We obtain a good cor- 
rection of the residuals, with a centroid at £0 = 6.27 keV and a 
FWHM o- = 0.07 keV. The flux of the Compton shoulder mod- 



elled in this way(5.8!^-^xlO""ph cm 
with that measured with S her pa. 



s ) is fully compatible 



4. Discussion 

4. 1. The nature of CG X- 1 

Off-nuclear Ultra-Luminous X-ray sources (ULXs) have been 
observed in several nearby galaxies by Einstein, ROSAT, 



ASCA and now by Chandra (e.g. [Colbert & Mushotzky 


1999; 


Mizuno et al. 


2001 


; pa Parola et al. 


2001 




Fabbiano et al. 



200 1| ). Their luminosities are in the range 10 - 10 erg 
s largely super-Eddington for stellar mass objects. Even if 
some of them can be expla ined as the product of young su- 



pernova remnants (e.g. [B01| , and references therein), the basic 
idea is that such luminosities are achieved through accretion 
onto black-holes with masses of 50-500 M0, i.e. intermediate 
between the masses of BH in Galactic binaries and those in 
AGNs. 

CG X-1 was first obse rved by ROSAT/HRl in 1995 (source 
3 inpuainazzi et al.||1999|), with a 0.5-2 keV flux of ^ 2 x 10 



Table 4. Count rates (0.5-2 keV) of CG X-1 for all 
ROS AT/HRl observations. The 0.1-10 keV fluxes are inferred 
from the powerlaw model parameters of the Chandra June 
2000 observation. 



Date 


Exp. (s) 


c/s 




Flux 0.1-10 keV (cgs) 


SEP95 


4193 


0.0045 


±0.0010 


9.5 X 10^1-' 


FEB96 


1081 


0.0019 


±0.0013 


4.0 X 10-13 


SEP96 


1850 


0.0011 


± 0.0008 


2.3 X lO-i-'' 


MAR97 


26673 


0.0023 


± 0.0003 


4.8 X 10"" 


AUG97 


46367 


0.0027 


± 0.0002 


5.7 X 10-" 



erg cm"^ s"', fully consistent with that of the low state of the 
June 2000 Chandra observation. The Circinus Galaxy was ob- 
served again by ROS AT/HRl in the following years: we down- 
loaded the data from the archive and analysed the count rates of 
CG X-1 for each of them (Table There is evidence of some 
variability, but the source seems always to be caught in a low 
flux state. On the other hand, a short-term, periodic variabil- 
ity is observed in the observations with the highest exposure 
times. The clearest evidence comes from the longest observa- 
tion, in August 1997, whose lightcurve power spectrum shows 
a peak at a frequency corresponding to a period of ^ 27 ks, 
as observed by BOl in Chandra data (see below). The 1990 



ROSAT/PSPC observation could not be used to study CG X-1 
because of its poorer spatial resolution (^ 25"), which is ac- 
tually larger than the distance between the source and the nu- 
cleus. 

CG X-1 was then ob s erved b y Chandra and s t udied i n de- 
tail by [Smith & Wilson| ( |200l| ), [Sambruna et~aL| ( |2001a[ ) and, 
in particular, BOl . Indeed, the best spectral information come 
from the longest observation, in June 2000, which also gives 
the cleares t evid ence of a short-term, periodical behaviour of 
the source (BOl ). Therefore, we will mainly refer to the results 
presented by ^01| and to their proposed interpretations of the 
nature of CG X- 1 , adding further insights using our analysis of 
all Chandra observations (see Table pt). Moreover, since it was 



shown in Sect. 3.2 that this source is likely to be the origin of 
most of the flux excess in the second BeppoSAX observation 
(with the exception of the 6.7 keV line and the PDS variation: 
see next section), we can assume that the spectrum resulting 
from the subtraction of the two BeppoSAX observations is a 
genuine spectrum of CG X-1. This spectrum can give valuable 
information, because, together with the Chandra June 2000 ob- 
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servation, is the only one which includes more than a period of 
the source. 

There are several pieces of evidence which strongly suggest 
that CG X-1 lies in the disk of the Circinus Galaxy (BOl). At 
the distance of 3.8 Mpc, the flux peak measured in March 2000 
corresponds to a luminosity of ~ 10"*° erg s this implies an 
object with a mass M > 80 Mq, if radiating at or below its 
Eddington limit. 

A n opt ical counterpart is not detectable down to my < 25.3 
mag ( BOl ). One of the most characteristic feature of CG X-1 
is its clear periodic behaviour: in the June 2000 observation it 
was found to display a flux variability by a factor ^ 20 in a 
period of 7.5 + 0.2 hr ( BOl ). However, long-ter m var iability is 
also observed, as already pointed out in Sect. 3.1.2: the flux 



peak in the June observation, for example, is about 1/4 of the 
observed mean in March. Periodic behaviour was also found in 
our analysis of the high flux state of March 2000, but the shorter 
exposure time prevented from confirming the above-mentioned 
period. The shape of the lightcurve during a period could be 
simply explained with an eclipse by a binary companion, but in 
order to make this model self-consistent an unphysically large 
companion would be required ( BOl ). 

is the mod- 



Another possibility, again suggested by BOl 



ulation of the accretion rate from accretion-disk instabilities, 
as observed in the quasi-periodic episodes of the Galactic mi- 



croquasar GRS 1915H-105 (3reiner et al. 1996; Belloni et al 



1997). The statistical quality of the available observations 



of CG X-1 is not good enough to perform a spectral vari- 
ability analy sis as detailed as the one performed by Belloni 
et al. (1997), which led them to explain the variability in 
GRS 1915H-105 in terms of cyclic al em ptying and refilling of 
the inner accretion disk. However, BOl claim a possible vari- 
ation of the spectral shape during the different phases of the 
periodical lightcurve and the different temperatures needed in 
our fits of all Chandra observations of CG X-1 (see Table |^ 
could be due to the fact that we are observing the source in 
different phases of the period. 

Regarding the long-term variability of CG X-1, if we com- 
pare the only two observations with a spectral shape averaged 
over more than a period (BeppoS AX and Chandra June 2000), 
we see that, despite a change in flux of a factor ^ 6, the temper- 
atures for the disk black body model are consistent with each 
other (see Table H). This suggests that the long term variability 
is not associated with significant variations of the inner radius 
of the accretion disc. 



of CG X-1, also when the source was as bright as supposedly 
was in BeppoSAX. 

We re-analysed data from the old BeppoSAX observation 
and added a line at 6.67 keV: the upper limit for its flux is 2.6 x 
10"^ ph cm^^s ' , which is indeed a factor of 3 lower than in the 



new data (Sect. 3.2. 1 ). On the other hand, the above-mentioned 



upper limit is fully consistent with the flux observed by the 



ph cm 



-2^-1. 



Sambruna 



gratings aboard Chandra (2.7 x 10" 
et al. 2001b) in the nuclear region. 

Contamination from another off-nuclear source is a possi- 
ble explanation. CG X-2 (see Fig. 1 in BOl) is an object with 
a flux comparable to that of CG X-1 in its low state, but fairly 
constant during all Chandra observations, being Fo.i-iokeV = 
1.3 X 10"'^ erg cm"^ s"'. Since it was completely absent in the 
1995 R OS AT IHKl observation and has optical and radio coun- 
terparts, BOl suggested it could be a young SNR. Our analysis 



of the other ROS AT IWSl observations shows the source is not 
present up to at least September 1997, i.e. six months before 
the first BeppoSAX observation. 

CG 



As observed by BOl 

(EW: 



X-2 has a huge iron line 
= 1.59 keV) at ^ 6.8 keV which is likely the blend of two 
lines (at 6.7 and 6.97 keV), as suggested by its intrinsic width of 
190 eV. The flux of the line, which we calculated in the longest 
Chandra observation, is 1.6 ± 0.4 x 10"^ ph cm"^ s"', which 
is almost a factor 5 lower than required in the new BeppoSAX 
data. The shorter exposure times of the other observations do 
not allow to search for variability of the flux of the line. 

Therefore, it is not possible to exclude that the ionized iron 
line observed in the second SAX observation, even if clearly 
including the Une of CG X-2 that was probably not present in 
the previous observation, could be an effect of a variation in the 
nuclear emission. 

The same explanation can be adopted for the variation in 
the PDS spectral shape. Although the total count rates are con- 
sistent with each other (see Table |l]), there is clearly a decrease 
in the lower part of the PDS band (Fig. ||). A possible inter- 
pretation is in terms of a change of the column density of the 
absorbing matter (the torus). However, it must be noted that the 
PDS field of view is very large (1.5x1.5 degrees) and so the 
emission from the nucleus could be contaminated by a num- 
ber of sources in an unpredictable way. Among these sources, 
an important contribution could come from the second source 
present in the LECS and MFCS images (see Sect. [3.1.1[ ): it may 
be that its flux over 10 keV was higher in the first SAX obser- 
vation. Unfortunately, there is no way to test this hypothesis. 



4.2. The ionized iron line and the PDS variation 

The spectral and flux properties of CG X- 1 are clearly the best 
explanation for the variations observed in the Circinus Galaxy. 
However, no iron line is present in the Chandra spectrum of 
CG X-1, while large residuals are clearly detected in the new 
BeppoSAX data in excess of the iron line observed in the old 
data. The centroid energy of the residuals is strongly shifted 
to 6.7 keV, indicating ionized iron. The flux required by the 
SAX data is much larger than the upper limit in the spectrum 



5. Conclusions 

The most likely explanation of at least most of the flux and 
spectral variation observed between two BeppoSAX observa- 
tions taken almost three years apart is in terms of a high-flux 
state of CG X-1, which is located well within the adopted 
BeppoSAX source extraction region of 2'. There are several 
pieces of evidence: 

- CG X- 1 shows a strong, long-term variation, reaching a flux 
which is comparable with that of the nucleus and is consis- 
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tent with that measured in the residuals of the two SAX 
observations. 

- The spectral shape of CG X-1, when averaged over more 
than a period, is fully consistent with the spectrum of the 
LEGS and MEGS BeppoSAX residuals. 

- Finally, the strongest evidence comes from the short-term 
variation of GG X-1, on a period of 27 ks: the MEGS 
lightcurve of the new observation clearly varies with the 
same period, indicating that a significant part of the ob- 
served flux originates from GG X-1. The periodic be- 
haviour was also present in the data from the old SAX ob- 
servation, but the variation amplitude was much less, be- 
ing consistent with contamination from a low-flux state of 
GG X-1. 

The best interpretation for the nature of GG X-1 is in terms 
of a > 80 Mq black hole in an accreting binary system in 
Gircinus, as previously suggested by BOl. Our analysis of all 



Chandra observations and the BeppoSAX spectrum shows that 
the long-term variability of this source is not associated to spec- 
tral variations, indicating that it is not due to changes of the 
inner radius of the disk. 

However, GG X-1 cannot be the cause of the variation of 
the PDS and the presence of an ionized iron line in the residu- 
als between the two BeppoSAX observations. In both cases, 
a variation of the properties of the circumnuclear matter in 
the AGN environment is a possible explanation: an increase 
of the column density of the torus, for example, would cause 
the observed decrease of the PDS count rates at lower energies. 
However, it is difficult to imagine the physical cause to support 
such an increase. 

As for the ionized iron line, it should be noted that the SAX 
residuals must include a line from GG X-2, which was likely 
not present in the old observation. Nevertheless, the contribu- 
tion of the flux of this line is too low to exclude a variation of 
the 6.7 keV iron line originating from the nucleus. 

Finally, our new analysis of the Chandra HEG nuclear 
spectrum has led to the detection of a Gompton shoulder in the 
6.4 keV iron line. Its flux is about 20% the line core, in agree- 
ment with theoretical expectations for Gompton-thick matter, 
providing one more argument in favour of the association of 
the matter producing the iron line with the 4 x 10^"* cm"^ neu- 
tral absorber. 
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